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Abstract. This article studies a hyperbolic conservation law that models a 
highly re-entrant manufacturing system as encountered in semi-conductor pro- 
duction. Characteristic features are the nonlocal character of the velocity and 
that the influx and outflux constitute the control and output signal, respec- 
tively. We prove the existence and uniqueness of solutions for L^-data, and 
study their regularity properties. We also prove the existence of optimal con- 
trols that minimizes in the L'^-sense the mismatch between the actual and a 
desired output signal. Finally, the time-optimal control for a step between 
equilibrium states is identified and proven to be optimal. 

1. Introduction and prior work. This article studies optimal control problems 
governed by the scalar hyperbolic conservation law 

dtpit,x) +da:iX{W{t))p{t,x)) ^0 where W{t) = [ p{t,x) dx, (1) 

Jq 

on a rectangular domain [0, T] x [0,1] or the semi-infinite strip [0, oo) x [0,1]. We 
assume that A(-) G C^([0, -|-cxd); (0, +oo)) in the whole paper. 

This work is motivated by problems arising in the control of semiconductor man- 
ufacturing systems which are characterized by their highly re-entrant character, see 
below for more details. In the manufacturing system the natural control input is 
the influx, which suggests the boundary conditions 
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p{0,x) = po{x), for < X < 1, and p{t,0)X{W{t)) = u{t), for t > 0. (2) 



Various different choices of the space of admissible controls are of both practical and 
mathematical interest, each leading to distinct mathematical problems. Motivated 
by this application from manufacturing systems, natural control objectives are to 
minimize the error signal that is the difference between a given demand forecast 
yd and the actual out-flux y(t) = X{W{t))p{t,l). An alternative to this problem 
modeling a perishable demand, is the similar problem that permits backlogs. In that 
case, the objective is to minimize in a suitable sense the size of the different error 



while keeping the state p{-,x) bounded. This article only considers the problem of 
perishable demand and the minimization in the L^-sense. 

Partial differential equations models for such manufacturing systems are moti- 
vated by the very high volume (number of parts manufactured per unit time) and 
the very large number of consecutive production steps which typically number in 
the many hundreds. They are popular due to their superior analytic properties and 
the availability of efficient numerical tools for simulation. For more detailed discus- 
sions see e.g. [5, 2, 3, 6, 16, 17, 19]. In many aspects these models are very similar to 
those of traffic flows, compare e.g. [12]. 

The study of hyperbolic conservation laws, and especially of control systems 
governed by such laws, have a rich history. A modern introduction to the subject 
is the text [8]. From a mathematical perspective, the choice of spaces in which 
to consider the conservations laws (and their data) provides for distinct levels of 
challenges. Fundamental are question of wellposedness, regularity properties of 
solutions, controllability, existence, uniqueness and regularity of optimal controls. 
Existence of solutions, regularity and well-posedness of nonlinear conservation laws 
have been widely studied under diverse sets of hypotheses, commonly in the context 
of vector values systems of conservation laws, see e.g. [4,7,9]. Further results 
on uniqueness may be found in [11], while [10] introduced an a distinct notion 
of differentiability of the solution of hyperbolic systems. For the controllability 
of linear hyperbolic systems, see, in particular, the important survey [22]. The 
attainable sets of nonlinear conservation laws are studied in [1, 15, 18,20,21], while 
[14] provides a comprehensive survey of controllability that also includes nonlinear 
conservation laws. 

This article is, in particular, motivated by the recent work [19] which, among 
others, considered the optimal control problem of minimizing \\y — yd\\L^{o,T) (the 

norm of the difference between a demand forecast and the actual outflux). That 
work derived necessary conditions and used these to numerically compute optimal 
controls corresponding to piccewise constant desired outputs yd- 

The organization of the following sections is as follows: First we rigorously prove 
the existence of weak solutions of the Cauchy problem for the conservation law (1) 
for the case when the initial data and boundary condition (2) lie in L^(0,1) and 
i^(0,T), respectively. Next we establish the existence and uniqueness of solutions 
for the optimal control problem of minimizing the L^-norm of the difference between 
any desired L^-demand forecast yd and actual outflux y{t) = X{W{t)) ■ p{t,l). 
Finally, in the classical special case where 



signal 




(3) 
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we prove that the natural candidate control for transferring the system from one 
equilibrium state to another one is indeed time-optimal. 

While preparing the final version of this article, the authors received a copy of 
the related manuscript [13] which is also motivated in part by [5, 2, 19] and which 
addresses wellposedness for systems of hyperbolic conservation laws with a nonlocal 
speed on all of K" . It also includes a study of the solutions with respect to the initial 
datum and a necessary condition for the optimality of integral functionals. There 
are substantial differences between [13] and our paper, especially the treatment of 
the boundary conditions and the method of proof. 

2. Existence, uniqueness, and regularity of solutions in L^. 

2.1. Technical preliminaries and notation. For any A G C^([0, +oo); (0, +oo)) 
define the functions A, A G C°([0, oo); (0, oo)) and d G C°([0, oo); [0, cx))) with re- 
spect to A as 

A(M) := inf X{W), X{M) := sup X{W), d{M) ~ sup \X'{W)\. (5) 

0<W'<A/ 0<H'<A/ 0<W<M 

For convenience we extend A to all of M in such a way that this extension, still 
denoted A, is in C^{R; (0, +oo)). 

2.2. Weak solutions of the Cauchy problem. First we recall, from [14, Section 
2.1], the usual definition of a weak solution to the Cauchy problem (1) and (2). 

Definition 2.1. Let T > 0, po e ^^(0,1) and u G L^iO,T) be given. A weak 
solution of the Cauchy problem (1) and (2) is a function p G C°([0, T]; L^(0, 1)) 
such that, for every r G [0, T] and every (p G (^""^([O, r] x [0, 1]) such that 

(^(r,a;) = 0,Va; G [0,1] and ^^(t, 1) = 0, Vt G [0, r], (6) 



one has 



pit, x) [ipt (i, x) + X{W{t))ip^ [t, x))dxdt 

u{t)ip{t,Q)dt + I po{x)ip{0,x)dx ^ 0. (7) 



"'0 



^0 

One has the following lemma, which will be useful to prove a uniqueness result 
for the Cauchy problem (1) and (2). 

Lemma 2.2. If p G C°([0, T]; L^(0, 1)) is a weak solution to the Cauchy problem 
(1) and (2), then for every t G [0,T] and every ip G C"'^([0,t] x [0, 1]) such that 

^(i,l) = 0,VtG [0,t], (8) 

one has 

1 pT 

p{t,x){ipt{t,x) + X{W{t))(p^{t,x))dxdt + / u{t)ip{t,0)dt 
Jo JO 

1 .1 

p{T,x)(p{T,x)dx + / po{x)ip{0,x)dx = 0. (9) 

JO 
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Proof. The case r = is trivial. For every t e (0, T] and e G (0,t), let rje G 
C^{[0,t]) be such that 

?7e(T) = and i]e{t) = 1, Vte [0,t - e] and ?7;,(t) < 0, £ [0,r]. (10) 

It is easy to prove that, for every h E C"([0, r]), 

lim / ri'^{t)h{t)dt = -h{T). (11) 

Then, for every e C^([0, r] x [0, 1]) satisfying (8), let (p^(t, x) := ri^{t)ip{t, x). This 
obviously verifies 

iPe{T,x) ^0,yx e[0,l] and ^^(t, 1) 0, Vt e [0, r]. (12) 

Since p € C°([0, T]; L^(0, 1)) is a weak solution to the Cauchy problem (1) and 
(2), we have 

rT /.I 

p{t,x){{^e)t{t,x) + X{W{t)){^,),{t,x))dxdt 



^0 



"'0 



+ / u{t){ip,){t,0)dt+ f po{x){ipe){0,x)dx = 0. (13) 

Jq Jo 
Using the definition of ip^, (10) and (13), one has 

p{t, x){ipt{t, x) + X{W{t))(pa,{t, x))dxdt 
/•I 

u{t)ip{t,0)dt + / pQ{x)ip{0,x)dx 
Jo 

1 

(1 - Ve{t))p{t, x){^t{t, x) + X{W{t))^^{t, x))dxdt 
[l ~ ri,{t))u{t)Lp{t,Q)dt - ( ( 7j'^{t)p{t,x)ip{t,x)dxdt. (14) 

t Jr-e Jo 

Observing that p G C°([0, T]; 2.1(0, 1)), A G Ci(]R; (0, cx))) and G Ci([0,r] x [0,1]), 
we point out that the functions W{-) ~ p{-, x)dx, p{-, x)ip{-, x)dx and \{W{-)) 
are ah in C°([0,r]). 

We can estimate the first two terms on the right hand side of (14) as 

r-l 

{\-rie{t))p{t,x){'pt{t,x) + \{W{t))if^^.{t,x))dxdt <Ke, (15) 

Jt-s Jo 
and 

{l-T]e{t))u{t)Lp{t,0)dt < K u{t)dt, (16) 

where if is a constant independent of e. While for the last term on the right hand 
side of (14); we get from (11) that 

r]'^{t)p{t, x)ip{t, x)dxdt = / '?e(i)(^ / p{t, x)ip(t, x)dx^ dt 



T-e JO 



p{t, x)ip{T, x)dx as £ — > 0. (17) 
In view of (15)-(17), letting e -> in (14) one gets (9). □ 
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Theorem 2.3. // po g L^(0,1) and u S L^(0,T) are nonnegative almost ev- 
erywhere, then the Cauchy problem (1) and (2) admits a unique weak solution 
p G C*'([0, r]; i^(0, 1)), which is also nonnegative almost everywhere in Q ~ [0, T] x 
[0,1]. 

Proof. Wc first prove the existence of weak solution for small time: there exists 
a small S G (0, T] such that the Cauchy problem (1) and (2) has a weak solution 
p G C°([0, (5]; i^(0, 1)). The idea is to find first the characteristic curve ^ = ^{t) 
passing through (0,0), then construct a solution to the Cauchy problem. 
Let 

^s.M {e G C"{[0,6]): ^(0) = 0,X{M) < ^^''^ ~ f^^ < A(M),Vs,t G [0,6] s> t}, 

(18) 

where A, A are defined by (5) and 

AI := |Kt||Li(o,T) + IIpoIIlmo,!)- (19) 

We point out here that the case d{M) = (by (5), A is a constant in [0, Af]) is 
trivial. We only prove Theorem 2.3 for the case d{M) > 0. 
We define a map F : ns,M C"([0, S]), ^ ^ F{(), as 

F{0{t) / A( / u{a)da + / po{x)dx)ds, G Qsm, V< G [0,6]. (20) 



Jo Jo Jo 
It is obvious that F maps into Qs,m itself if 

< (5 < T and J < =J—. (21) 
A(M) ^ ' 

Now we prove that, if d is small enough, is a contraction mapping on fls,M with 
respect to the C" norm defined by 

||Cllc°([o,5]) := sup [^{t)[. 

0<t<S 

Let a, 6 e ^SM-We define I, G C°{[0,S]) and G C°([0,(5]) by ^t) := 

max{a(t), a(t)} and ^2i^) •= niiii{a(i), ^(t)}- By (5) and changing the order of 
the integrations (see Figure 1), we have 

[F{^2m - F{^,){t)[ < d{M) / / po{x)dx ds 

Jo Jl-Ci(s) 

= d(A/) / _ po{x)it-^^ {l~x))dx 

+d(M) / _ po(a:)(?^'(l-x)-$7\l-x))dx. 

< d{M) f p^{x)dx ■ i^^'iUt)) - ii\um 

+d{M) _ po{x)(^2\^ - ^) ^Ii'i'^ - x))dx 

< d{M) po{x)dx- sup {l2\y) -Ii\y))- (22) 
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Using the definitions of ^1,^2 ^-nd of fls,M, wc obtain that, for every y G [0, ^2(0] 
(sec Figure 2), 

- (?2 - ;^^"(-^) ) + '(^) ^-Ciy) 
< ^f,-e2(^^^^^^^^)U^fci(^^^^^^^^' 

- XiM)^ 2 'J A(Af)V^'^ 2 
1 /7 .?^'(^/)+?2"'(^/)^ 7 .?^'(y)+?2"'(2/)^ 



A(Af) 
Therefore, 



A(A/) 

< J—\\^, - ^.Wco^^oM)- (23) 



< 



A(M) Ji-c,(t) 



A(M) 



po(a;)dx ■ ||Ci - 6llc«([o,5])- (24) 



5; SjCO) 




1-5,(1) l-5j(t) 




|;| C'(y)+C'(y) j y ^| ^:'(y)+^:'(y) j 5,(1) y 



Figure 1. Change 
order of integrations 
for X and s in (22) 



Figure 2. Using 
property of fls m in 
(23) 



Since po £ -^^^(0, 1), we can choose S £ (0, 1) smah enough such that 

Po(x) ax < 



2d{M) ' 



(25) 



Then 
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- i^(6)llco([o,^]) < - e2||co([o.^])- (26) 

By means of the contraction mapping principle, there exists a unique fixed point 
^ = F(^) in ^ls,M- By (20), the fix point ^ is an increasing function in C^([0,(5]), 
and one has 

ft fi-m 
e'W = A(/ u{a)da+ Po{x)dx), yte[0,S]. (27) 

Jo Jo 

Then we define a function p by 

(po{x-m), 0<e(0 <2:< l,0<t <(5, 

which is obviously nonnegative almost everywhere. Direct computations give that, 
for every t e [0,5], 

0<W{t) = [ p{t,x)dx 







dx + po{x — ^{t))dx 



< 



ei^-'im-^)) Jm' 

u{a)da+ / Po{v)dy 
Jo 

\\u\\Li(^o,T) + \\po\\moA)^M. (29) 



Using (5), (27) and (29), we obtain the following estimates of ^' from above and 
below: 

< A(A/) < C'(t) = X{W{t)) < X{M), Vt G [0,5]. (30) 
We now prove that p e C°([G, S];L'^{0, 1)). For every s, t e [0, S] with s >t, 
1 

\p{s, x) — p(t, x)\dx 







< 



eic-'m-^)) e{^-'im'^))\ 

+ [pis,x) - pit,x)[dx+ [poix - ^{s)) - po{x ~ ^{t))\dx. (31) 

Jm Jm) 

As for the first term on the right hand side of (31), we choose C 
C^([0, T]) which converges to u in L^{0, T), then we have 
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?(*) 



< 



dx 



< 



+ j )|u(cr) -u"(cr)|dcr 



dx 

dx 
dx 



< 2 



\u{(j) - w"(cr)|(icr 



dx 



dx. 



By (30), 



m 



< 



dx 



«W w"(^-i(^(s)-x))-u"(r^(e(i)-2:)) 



1 



u(a)dadx 



+c„ 



Poiy)dydx, 



(32) 



dx 



(33) 



where C„ is a constant independent of s and t but depending on u". By changing 
the order of integrations, we obtain furthermore (see Figure 3) 



u{a)dadx 



< 



' u{a)da-\as)-m\ 







< MLHo.T)-\m~m\ 



u{a)dxd<7 



(34) 



and (see Figure 4) 
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/ / Po{y)dydx 

Jo Jl-i{s)+x 

/ / + / / + / / )poiy)dxdy 

< ( / + / + / )poiy)dym-m\ 

< \\po\\LHos)-\as)-m\- (35) 




As for the second term on the right hand side of (31), it is easy to get that 



\p{s,x) — p{t,x)\dx < / p{s,x)dx+ / p{t,x)dx 

u{(7)da+ / Po{y)dy. (36) 

^0 



As for the last term on the right hand side of (31), we choose {Po}^=i C C^([0, 1]) 
which converges to po in L^{0, 1), then we have 
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1-1 

\po{x - - po{x - £,{t))\dx 

as) 

< [ \poix ~ a--^)) - p^ix - ^is))\dx + [ \p^{x-as))-p'Six-m)\dx 

+ [ \p'six^m)-poix-m)\dx 

/.l-?(s) ri-m . 

< i + )\po{y)-p'i{y)\dy + D„\as)-m\ 

< 2/ \po{y) - pUy)\dy + Dn\m ~ mi (37) 

Jo 

where D„ is a constant independent of s and t but depending on Pq. 

Using (19) together with the estimates (31) to (37), we obtain for any s, < G [0, 6] 
with s >t, 

1 

p(s, x) — p{t, x)\dx 

< 2/ \uia)-u^ia)\da + C„|e(s) - CWI + A/|^(s) - S(i)l 
Jo 

r'(«M-?(t)) rm-m 

u{a)da+ / Po{y)dy 
Jo 

+2 [ \p,{y)-p^(y)\dy + D„\as)-m\- (38) 
Jo 

We can choose u" and such that |w(cr) — u"(cr)|(icr and \po{y) ^ Po{y)\dy 
are small as we want. Then according to (30) and the fact that u G L^{0,T) and 
po S £""^(0,1), the right hand side of (38) is sufficiently small if s and t are close 
enough to each other. This proves that the function p defined by (28) belongs to 
CO([0,5];LH0,l)). 

Next, we prove that p defined by (28) is a weak solution to the Cauchy problem 
(1) and (2). Let 6' e [0,6]. For any ip G C^{[0,S'] x [0,1]) with (p{d',x) = and 
ip{t, 1) = 0, let 

/ / pit,x){iptit,x) + \{Wit))ip^{t,x))dxdt. (39) 



Then we have 

5' nl 

/ po{x - ^{t)){^t{t,x) + X{W{t))^,{t,x))dxdt 

Jm 

S' rt 

u{<j){^t{t,m - m)) + x{w{t))^At,m - m))d<jdt 



Jo 



s' i-i-m 

Po{y){^t{t,m + y) + HW{t))^At,m + y))dydt, (40) 



JO 



and thus 
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u{a)ipia,0)dt - [ poivMO, y)dy. (41) 
Jo 

This proves the existence of weak solutions to the Cauchy problem (1) and (2) for 
small time. 

Now we turn to prove the uniqueness of the weak solution. Let us assume that 
p € C^{[0, S]; L^{0, 1)) is a weak solution to the Cauchy problem (1) and (2). Then 
by Lemma 2.2, for any t G [0,(5] and ip e C1([0,t] x [0,1]) with -ipit, 1) = 0, 

p{t,x){Mt,x) + HW{t))ipx{t,x))dxdt + [ u{t)^P{t,0)dt 



^0 -'0 

- / p{T,x)ip{T,x)dx + / po{x)ip{0,x)dx = 0, (42) 
^0 Jo 

where W{t) := J^p{t,x)dx. 

Let ^(i) := X{W{s))ds and ipo S Cg(0,l) (i.e. a function with compact 
support in (0, 1)). Then we choose the test function 

>^ fV'0(?(T)-CW+^), 0<X<l{t)-l{T) + l,0<t<T, 

\0, 0<^{t)-^iT) + l<X<l,0<t<T, ^ ' 

which obviously belongs to C^([0,r] x [0,1]) and satisfies the following backward 
Cauchy problem: 

>t + A(T^(i))^:, = 0, 0<t<T,0<.T<l, 

i}{T,x)=i!o{x), 0<a:<l, (44) 

1)=0, 0<t<T. 

In view of (42), we compute 
1 

p(t, x)-iljQ{x)dx 

.l-C(r) 

u{t)ijo{^{T) - i{t))dt + / Po{x)Mi{r) + x)dx 

Jo 

^=M^^^^My)dy+ f po{y~l{T))My)dy. (45) 
^0 a? (e(r)-y)) hir) 

Since V'o S ^^0(0, 1) and t G [0, S\ were arbitrary, we obtain in C°([0, J]; L^(0, 1)) 
that 

'p^{x-l{t)), 0<l{t)<x<l,Q<t<5, 



which hence gives 
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= M p{s,x)dx)ds 
Jo Jo 

A( / u{a)da+ / po{y)dy)ds 
Jo Jo 

= nm)- (47) 

It is easy to check that ^ e ^s.m when 5 is smah enough, which impUes that ^ = ^ 
since ^ is the unique fixed point of F in fls^Ai for 6 small enough, and then p = p 
by comparing (28) and (46). This gives us the uniqueness of the weak solution for 
small time. 

Now we suppose that we have solved the Cauchy problem (1) and (2) to the 
moment t G (0,T) with the weak solution p e C°([0, r]; Li(0, 1)). By (29), the 
following uniform a priori estimate holds for every t £ [0, r]: 

< W{t) = I p{t, x)dx < M. (48) 
Jo 

Hence we can choose 5o e (0,T) small enough such that (21) holds and 

,_,/(^'^)^^^MM)- 

Applying the previous results on the weak solution for small time, the weak 
solution p G C°([0, r]; L^(0, 1)) is extended to the time interval [r, t + Sq] Ci [r, T]. 

Moreover, by (48) and po S ^^(0, 1), u £ L^{0,T), one can find a suitably small 
5q > independent of t such that 

' A(AJ) ' " 



p{T,x)dx < sup / u{a)da + sup / po{y)dy 

i-So teloT — ^a-i-'t xe[a,i-5o]Jx 



X(M) , , 

^ im- ^''^ 

Step by step, we finally have a unique global weak solution p € C'^{[Q,T]; L^{0,1)). 
This concludes the proof of Theorem 2.3. □ 

2.3. Remarks. 

Remark 2.4. Let p be the weak solution in Theorem 2.3. Let W G C°([0,T]) be 
defined by W(t) := p(t, x)dx and let ^ G C^i[0, T]) be defined by requiring 

e(o) = o, m^x{w{t)),vte[o,T]. 

Then, it follows from our proof of Theorem 2.3 that 

f poix - e(i)), < at) <x<i,o<t< r Hi), 



p{t,x) = < 



U'(?-MCW-^))' 

Moreover, W{t) can be expressed as 
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1 \ I u{a)da+ Poiy)dy, < t < ^'(1), 

W{t)= I p{t,x)dx= {■'jt ■'^ (51) 



which inipHcs that 



< W{t) = [ p{t, x)dx < M, \/t G [0, T] 
Jo 



(52) 



and 

< A(M) < e'(i) = X{W{t)) < \(M), yt e [0, T]. (53) 
Finally, W is absolutely continuous: 

W{t) = W{0) + [ W'{s)ds (54) 







with 



(u{t)-eit)poii-m), o<f <r'(i), 
V^'^' e(5-(^(i)-i)) ' '-^ 

and 

< / \W'it)\dt< M. (56) 

Remark 2.5. (Hidden regularity.) From the definition of the weak solution, 
we can expect p e i^(0, l;L^(0,r)). However, the weak solution is more regular 
than expected. In fact, under the assumptions of Theorem 2.3, we have the hidden 
regularity that p e C"([0, 1]; ^^(O, T)) so that the function t ^ p{t,x) G L^{0,T) 
is well defined for any fixed x G [0, 1]. The proof of the hidden regularity is quite 
similar to our proof of p G C°([0, 1]; L^(0, T)) by means of the exphcit expression 
of p (sec also (54)-(56) that we use when T is large). 

Remark 2.6. If po G iP(0, 1) and u G LP{0,T) {p > 1) are nonnegative al- 
most everywhere, then the Cauchy problem (1) and (2) admits a unique weak 
solution p G C°([0,T];LP(0,1)) n C°([0, 1]; iP(0, T)), which is also nonnegative 
almost everywhere in Q = [0,2^] x [0,1]. In fact, the uniqueness of the weak 
solution comes directly from Theorem 2.3. And the expression of the solution 
pG CO([0,r];ii(0, l))nC"'([0,l];ii(0,T)) given by (50) shows that p belongs to 
CO([0,T];LP(0,l))nC"([0,l];iP(0,r)). 

Remark 2.7. If po G C^([0, 1]) and u G C^([0, T]) arc nonnegative with 
'«(0)-po(0)-0 

(57) 

u'iO) + X{ poix)dx)p'a{0)^0, 



then the Cauchy problem (1) and (2) admits a unique classical solution p G ([0, T] x 
[0, 1]), which is also nonnegative. 
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3. L^-optimal control for demand tracking problem. Let po £ ^^(0,1) be 
nonnegative almost everywhere and let T > be given. Let us define 

L^(0,T) := {u e L^{0,T);u is nonnegative almost everywhere}. 

According to Remark 2.6, for every u G i^(0,r), the Cauchy problem (1) and (2) 
admits a unique solution p £ C°{[0,T], L^{0,1)) n C°{[0,1], L^{0,T)). 

For any fixed demand signal yd S i^(0, T) and initial data pQ, define a functional 
onL^(0,T) by 

J{u):^ ( \u{t)\^dt+ [ \y{t)-yd{t)\^dt,ueLl{0,T), (58) 
Jo Jo 

where 

y{t) := p{t,l)X{W{t)) (59) 
is the out-flux corresponding to the in-flux u £ L^(0,T) and initial data pQ. 

Theorem 3.1. The infimum of the functional J in L'^(0,T) is achieved, i.e., there 
4( 



exists Uoo G L^(0, T) such that 



J(u^) = inf J(u). (60) 

ueLl{0,T) 

Proof. Let {wnj^^Li C L^{0,T) be a minimizing sequence of the functional J, i.e. 

lim J{un) ~ inf J{u). (61) 

Then we have 

\\Un\\mo,T) + \\yn\\mO,T) <C, Vn e Z+. (62) 

In (62) and hereafter, we denote by C various constants which do not depend on n. 

The uniform boundedness of u„ in (0, T) shows that there exists Uoo G i+ (0, T) 
and a subsequence of {unk}'k^=i such that u„^. ^ Uqo in L^{0,T). For simplicity, 
we still denote the subsequence as 

Let pn be the weak solution to the Cauchy problem of equation (1) with the 
initial and boundary conditions 



p(t,o)A(iy(t)) = u„(i), o<t<r, 

p(0,x) =po(a;), 0<a;<l. 



(63) 



Let Wn ■■ [0, T]^R and ^„ : [0, T] E be defined by 

Wn{t):= f Pn{t,x)dx, Ut)-= [ X{Wn{s))ds. (64) 
Jo Jo 

Thus by (51), we have 
Ut) ^ I kI Un{<7)da + I ^"^ ^ pQ{x)dx)ds, 0<t< min{e-'(l),r}. (65) 



In view of (62) and (64), we can derive from (51) that 

||W„||co([o,T]) <C, V71GZ+, (66) 
which in turn gives with (64) that 

||^„||ci([o,T]) < C, VnGZ+. (67) 
Moreover, let us point out that is uniformly bounded from above and below: 
< A(C) < Cit) = KWnit)) < X(C), yt G [0,T], Vn G Z+, (68) 
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where A, A arc defined by (5) with 

C := sup ||w„||li(o.t) + ilPollLi(o.i) < oo. (69) 

nez+ 

Then it follows from Arzela-Ascoli Theorem that there exists G C''([0,T]) 
and a subsequence {^„, such that — s- in C*'([0,T]). Now we choose the 
corresponding subsequence {ttmli^i and again, denote it as {un}^=i- Thus we have 

Un ^ Uoo in L^(0,T), as n — > oo (70) 

and 

in-^U in C"([0,r]), asn^K). (71) 

Then one has 

= \x- urJm = iCooie^o^^;)) - urJm ^ o, as « ^ ^ (72) 

uniformly for x G [O,|'oc(r)). Thus we get for any xa e [0, ^^^(r)), 

C'^?;.' inC"([0,xo]), asn^oo, (73) 

and therefore, by passing the limit n ^ oo in (65), 

_ rt /-i-Cools) ^ 

CooW = M u^{a)da + / po{x)dx)ds, < t < min{^^ (1),^}. 

Jo Jo Jo 

(74) 

Let Poo be the weak solution to the Cauchy problem of equation (1) with the 
initial and boundary conditions 

ipit,0)XiW{t))=u^{t), 0<t<T, ^^^^ 
\p{0, x) = Po{x), < a: < 1, 

and denote 

M^oo(<):= / Poo{t,x)dx, eoo(i):= / X{Woo{s))ds. (76) 
Jo Jo 

We claim that = Coo- In fact, 

eoo(i) = X{ Uooi<j)da + / poix)dx)ds, 0<t< min{C'(l), T}- 

Jo Jo Jo 

(77) 

As in the proof of Theorem 2.3, there exists S > small enough which is depending 
only on Uqo and po such that 

C(t) = ^co(e)(0 := / X{[ u^{a)da+ [ \^{x)dx)ds (78) 
Jo Jo Jo 

has a unique fixed point in fi^^ (replacing M by C in (18)). This implies from (74) 
and (77) that ^oo(i) = ^oo(^) on [0,5]. Moreover, with the help of (49), there exists 
(5o > independent of r G (0,T) such that if Coo(t) = S.ooi''') ^i^cn £,oo{t) = Coo(^) 
on [t,t + 5o] n [t,T]. 
Therefore 

eoo=?oo and Cn^^oo inC°([0,r]), as n -> (79) 
and it follows that 
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Wn^W^, C^^L inC°([0,T]), asn^oo (80) 

and, for any xq £ [0,^oo(r)), 

C'^C' inCO([0,xo]), asn~>«). (81) 

Next prove that y„(t) = X{Wn{t))p„{t, 1) converges to y^{t) = X{W^{t))p^{t, 1) 
weakly in 2.2(0, T). By (62), {y„}^=i is bounded in L'^{0,T). Hence, it is suffices 
to prove that for any g G C^([0, T]), 



lim / ivnit) ~ yooit))g{t)dt ^ 0. 



(82) 



If Coo(r) < 1, then C«(r) < 1 for n large enough. By (50), (64) and (76), for 
every Xq S [0,^oo(7")), we have 

{Vn{t) ~ yoo{t))g{t)dt 



(po(l - in{t))\{Wn{t)) - Po(l - S.oo{t))X{W^{t)))g{t)dt 



< 



(T) 

1 



po{y)g{C\^ ~ y))dy 



Po{y)g{^-Hi - y))dy 



l-«oc(T) 



l-?o=(T) 
l-«=c(T) 



Po(2/)(5(C'(l-2/))-<?(C'(l-y)))^^2/ 



po{y)9{in^{^-y))dy 



< C sup |Ci(x)-C-^(a;)|+C|^oo(T)-xo|Uc|^„(T)-eoo(r)|i (83) 

0<-X<.Xq 

By (79) and (81), it is easy to get (82) from (83). 

If eoo(T) = 1 (i.e., T = e-i(l) ), for every r e [0,^-i(l)), we have 



{ynit)-yocit))g{t)dt 



< 



{yn{i) - yoc{t))g{t)dt 

{yn{t) ~ yoo{t))g{t)dt 

J T 

{yn{t) - y^it))g{t)dt\ + C(C'(1) - T)i 

Since it is known that for every r £ [0, T) 

{yn{t) - yoo{i))9{t)dt ^0 as n ^ oo, 

one has (82) for T = ^ooHl) from (84). 

If S,oo{T) > 1, then ^n(T) > 1 for n large enough and we have 



(84) 



{yn{t) ~ yoo{t))g{t)dt = 



Hi) 



+ 



(ynit) ~ yooit))git)dt. (85) 



From the above study, we need only to estimate the last term in (85). Assuming 
Cn^(l) < 'fTO^(f) (the case Cn^(f) ^ ^oo^(I) can be treated similarly), we get from 
(50) that 
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iVnit) - yoo{t))g{t)dt 



iuniC\Ut) - 1)) - Uoo{C^\U{t) - mgit)dt 



^-(^^ u^{a)girj^iamL{a) 





da 
da 



-da 



da 



C^»('7oo(cr)) 
where we denote 

Tnit) ^-'iUt) - 1), '7«W inHUt) + 1), 
Too(t) := ^^\^oo{t) - 1), ryoo(t) := + 1). 

From (79) and (81), we get 

~* Too ill C''^([Oi 2^1)5 as n — > 00 
and, for every to e [0, roo(T)), 

Vn Voo in C°{[t, Oo]), as n ^ 00. 
Therefore, by (86), one has for every to e [0,Too(r)) and for n large enough, 



iVnit) - yooit))git)dt 



< 



Wn(CT)g(?7„(cr))C;(cr) 



da 



T^(T) JO ' ^n{^n{(y)) 

^"^"^^ ' Un{a)g{j]n{(yy)^n{^) Uoo{cr)g{Vooi(j))S.'ooi'^] 



10 ^ ^'niVnicr)) CLiVocicr)) 

< C|r„(T)-ro.(r)|Uc|r„(C(l))|^ 



da 



-{giVnicr)) - g{r]oo{cr))da 



(cr)) 

^n(^).9(^oo(a))(-f^ - jrf^, 

[Un[a) ~ Uoa[a)) — — — da 



da 



10 CooiVooicr)) 
< C\r„iT) -r^{T)\i+ C\T„i^^\l))\i + C\r^{T) - to|^ 



+C sup |?7„ (cr) - ?/oo(o-) I + C* sup 

0<a<to 0<<T<to 



^U'yr.H) ^^(^00 (a)) 



4^X3 Woo (ct)) 
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By (70), (80), (89)-(90) and the arbitrariness of io G [0,Toc(T)), we have (82) for the 
case £,oo{T) > 1. This concludes the proof of (82). 
As a result, 



T pT 



J(moo) = / \uoo{t)\'dt + J ly^t) ~ydit)\'dt 

< liminf / |u„(t)pdt + liminf / \yn{t) ^ yd{t)\^dt 



n — ^oo /n n — ^oo 







< liminf J(u„)= lim J(un) ~ inf Jiu) (92) 

Ji— foo ri^oo ueL^{0,T) 

This shows Uqc is a minimizer of J(u) in L'^{0, T), and it proves also that u„ tends 
to Uoo strongly in i^(0, T). □ 

4. Time-optimal transition between equilibria. In this section, we focus on 
the specific model that relates the nonlocal speed to the total mass according to 
the assumption (4). 

It is immediate that constant boundary data p(-,0) = pin > eventually drive 
the state to the equilibrium p = p^- Together with the symmetry {t, p{t, x)) — > 
{T — t,l — X, p{T — t,l — x)) of the conservation law (1) this establishes (long-time 
state) controllability. Of particular interest is the question of how long it takes 
to drive the system from one equilibrium state po to another equilibrium state pi, 
compare also the numerical studies of transfers between equilibria in [19]. 

We first explicitly calculate all quantities for the corresponding piecewise constant 
boundary data p(-, 0), and subsequently prove that this boundary control is indeed 
time-optimal. 

Suppose pi > Po ^ arc constant, the initial density is the equilibrium p(0, x) = 
Po for x G (0, 1], and the desired terminal density is p{T,x) = pi for x e [0, 1] and 
some minimal T > 0. The case po > Pi > is similar. 

A natural choice for the boundary values is p(t, 0) ~ pi for t > 0. This determines 
for < t < T the control influx and the outflux via u{t) = piX{W{t)) and y{t) = 
poA(M^(t)), where is a solution of the initial value problem 

This can be integrated in closed form, yielding 



W{t) = -l + v/(l+Po)2 + 2t(pi-po) (94) 

and similar expressions for the fluxes and the speed. All characteristic curves are 
translations of the solution of the initial value problem 

e{t) = xiw{t))= ] m = o. (95) 

which has the explicit solution 



^ v/(l + Po)^ + 2t(pi-po)-(l + po) ^gg^ 
pi - po 

The time T to achieve this transition between equilibria is uniquely determined by 

£,{T) = 1 and evaluates to 

T=l+P^±Pl. (97) 
In the sequel we prove that this time is indeed minimal. 
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Note that W is a. continuous function, and, in particular W{T) = pi. It is 
convenient to extend p, u, v, and W to negative times by setting p{t, x) — W{t) = po 
and u{t) ~ y{t) ~ i+pQ ^ "0 = 2^0 for all i < 0. Then u is continuous except for a 
jump at t = 0, and y is continuous except for a jump at T. Note that the height of 
the jump of u at < = is larger than the corresponding jump of y at T. 



u{0+)-u{0-) = 



Pi - Po 
1 



whereas y{T+) - y(T') = 



Pi - Po 
1 



(98) 



Po 1 + Pi 

Supposing a jump of the reference demand from yd{t) ~ poX{po) for i < T to 
yd{t) = piX{pi) ior t > T at this earliest feasible time, the total backlog at any 
t > T is, due to the inverse response, 



m 



(2/d(s) - y{s))ds 



PoT 



1 



Po 



po I X{W{s))ds 

10 

V ' 

=1 



poT - po- Pq 
l + Po 



(99) 



Using the expression (97) for T, this simplifies for t > T to 

{Pi - Po)po 



m 



1 + Po 



yo{pi - Po) 



(100) 



W2 






u{t) 

^^^^ U2 = J/2 




J 






U 1 





Figure 5. Total mass, influx, and outflux during optimal transi- 
tion between equilibria 



Correspondingly, for < i < T the total mass W{t) < pi continues to grow, and 
hence the speed is further decreasing. Therefore, the influx u{t) = piX{W{t))) is 
larger than the eventual new equilibrium influx ui ~ piA(pi). The total excess in 
influx evaluates to 

a{T) ^ f{piX{W{s)) piX{pi))ds = ^P^ - P^'^P^ ^ u,{p, - po) (101) 
Jo 1 + Pi 

Together with the nominal difference (yi — yo)T between the accumulated equilib- 
rium fluxes over the time interval [0,T], these add up the difference in total mass, 
compare the three shaded regions in Figure 5, 
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W{T)-W{0)= f {u{s)-y{s))ds^a{T) + p{T) + {p^X[p^)-po\{po))T = pi-po- 
Jo 

(102) 

While it may seem intuitive that this control is time-optimal, we need to rigor- 
ously prove that it is indeed not possible to improve on this time by e.g. temporarily 
increasing the speed via smaller influxes. 

Proposition 4.1. The minimum time to transfer the state from one equilibrium 
p(0,x) = Po, X G (0, 1] to the equilibrium p{x,T) ~ pi > pQ, x G [0, 1] using influx 
w G Li([0,cx)),[0,cx))) isT=l + ^^^. 



t 




Proof. Suppose T > and p{t, 0) is an integrable function on [0, T] such that the 
solution of (1) satisfies p{T, •) = pi. 

Since p is constant along the characteristic curves, there exists to G [0, T] such 
that for all t G [to,T], p{t,0) = pi. Let ^: [0, T] i-^ [0, oo) be the unique function 
satisfying ^(0) = and ^'{t) = Mio p{t,x)dx). Then there exists a unique ti G 
(0,r] such that ^{ti) = 1. 

For < t < fo, C'(^) is bounded above by 

^'(t) = < -; 7^ — (103) 

1 + J^pit,x)dx 1 + J^^^^p{t,x)dx l + Po(l-e(i)) 

Rewrite as ((1 + po) — Po(,{t))C{t) ^ 1 a-nd integrate from t ^ to t ^ to to 
obtain a lower bound for tg. 

(l + Po)e(to)-^M(io)' <to. (104) 
The primary interest is the case of to < ti. For to < t < ti estimate 
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1 + Pit, x)dx + J^^^^ p{t, x)dx 

1 



(105) 



i + pi(ew-e(to)) + Po(i-e(i))' 

and integrate from to to ti to obtain 

(1 + Po - Pom^mih) - iih)) + \{Pi~ Po)(C(ii) - ato)f < <i - h. (106) 
Analogously, for h <t<T, the bound C{t) < 1/(1 + Pi - C(io))) yields 

^T) - e(ti) + Ipi mT) atojf m) ato)f) <T-h. (lo?) 

After combining the estimates (104), (106), and (107), elementary simplifications 
yield 

T > (l + Po)e(io)-^PoC(to)' 

+(1 +P0- po^tomiti) - ato)) + lipi- pomti) - ato))^ 
+e(T) - c(ti) + ((e(r) - ato)f m) ato))') ■ (108) 

Noting that ^(ti) = ^{T) - ^{to) 1, (108) simplifies to 

T>l + ^^i^+e(to). (109) 
This shows that the optimal choice is = 0, i.e. p{t, 0) = pi for all t > 0. 

It remains to dispose of the case when ti < to- For < t < ti, use £,'{t) < 
1/(1 + po{l - ^(t))) and S,{ti) = 1 to obtain 

h > e(ti)(l + Po) - \po^{tif = 1 + Ipo- (110) 

Similarly, for to<t<T, use ^'(t) < 1/(1 + pi{m - ^(io))) and C{T) - ({to) = 1 
to obtain 

T~to> (e(T) - e(to)) + ^Pi(e(T) - e(to))' = 1 + ^Pi. (Ill) 
Combining (110) and (111) together with to > ti yields 

T = (T - <o) + [to - h) +h>{l + ipi) + + (1 + ipo) > 2 + P^±Pl. (112) 

This shows that any controls for which ti < to will perform even worse than the 
ones in the first case. □ 
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